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We report on a study of vortex pinning in nanoscale antidot defect arrays in the context of the 
London Theory. Using a wire network model, we discretize the array with a fine mesh, thereby 
providing a detailed treatment of pinning phenomena. The use of a fine grid has enabled us to 
examine both circular and elongated defects, patterned in the form of a rhombus. The latter display 
pinning characteristics superior to circular defects constructed with the similar area. We calculate 
pinning potentials for defects containing zero and single quanta, and we obtain a pinning phase 
diagram for the second matching field, H — 2$ . 
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PACS numbers: 74.20.-z, 74.60. Ge, 74.76, 
I. INTRODUCTION 

It has long been known that defects in superconduc- 
tors can act as pinning centers for vortices. The trap- 
ping of vortices reduces dissipation and permits higher 
critical currents, which is important in engineering ap- 
plications involving superconductors. In nature, one ex- 
pects defects to be randomly distributed within a sam- 
ple. To examine vortex pinning in a more controlled set- 
ting, a body of experimental work has examined pinning 
phenomena in artificial defect arrays. Advances in sub- 
micron lithographyEnl permit the construction of artifi- 
cial nanoscale defect arrays. Lorentz microscopy permits 
direct imagingEl of the magnetic fields associated with the 
vortices and supercurrents. Scanning Hall probesBE2l pro- 
vide less detail, albeit over a greater spatial range. The 
defects can be fashioned as regions of excised supercon- 
ductor, "antidots" , or they may consist of ferromagnetic 
inclusions. This study concentrates, on the antidot case. 

The recent experimental worktTB on nanoscale defect 
arrays in. thin films has motivated a number of theoretical 
studiesfLd~li-3 of pinning phenomena in them. In this vein, 
we offer a detailed numerical treatment of thin-film pq 
odic antidot arrays in the context of the London theory 
Our study differs from most previous theoretical work in 
that it explicitly treats the thin-film geometry, rather 
than a bulk system with two dimensional ordering. Al- 
though we report on results obtained for the triangular 
antidot lattice, our method can readily be adapted to 
other geometries such as the square lattice. It is also 
within the scope of our method to consider more general 
inclusions, such as nanoscale defects with a permanent 
magnetic moment. 

Since the London Theory does not admit a solution 
in closed form for our system, a numerical approach is 
needed. In our case, the numerical treatment involves 
solving the London equations in discrete form. In doing 
so, we model the supercurrents in our system as a wire 
network. To satisfy the zero-current constraint within 
the inclusions, we expand the current in terms of ele- 
ments of a special basis, to be described in greater detail 
in the next section. To determine the optimum current 
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configuration, we minimize the discretized version of the 
London free energy with respect to the basis expansion 
coefficients. 



Synopsis of the results 

As shall be discussed, our method permits us to com- 
pute the free energy for systems containing a fixed num- 
ber of vortices at predetermined positions. By computing 
the London energy for different vortex configurations, we 
calculate the pinning potential for the vortices. Pinning 
potentials computed for different antidot radii indicate 
that pinning centers with radii larger than £(T) are more 
efficient at pinning than those whose size is in the vicinity 
of the coherence length. This result is in agreement-with 
one of the experiments! as well as theoretical workHl on 
three dimensional superconductors containing columnar 
defects. The pinning potentials also allow one to identify 
interstitial pinning. With the aid of pinning potentials, 
we have produced a phase plot displaying the optimum 
number of pinned vortices as a function of defect radius. 
A salient feature is the weak dependence of the number 
of pinned quanta on the effective magnetic penetration 
depth, X e ff- Finally, we have examined the pinning po- 
tentials of defects with an elongated shape. To consider 
defects with sharp corners, we have studied antidots in 
the shape of a rhombus. Two notable features of the 
pinning potentials of rhombic defects are a lower pinning 
energy and a higher escape barrier than that of circular 
defects of comparable area, qualities indicating that the 
rhombus is superior to the circle as a pinning center. 

The article is organized as follows: In Section II, we 
present the wire network model for the supercurrents in 
the nanoscale defect array. Section III outlines our meth- 
ods for calculating the London free energies and current 
configurations, and Section IV reports our primary re- 
sults. We conclude in Section V with a summary. 
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II. THE MODEL 



To proceed, we consider antidots of radius R forming 
a regular triangular lattice in a thin-film superconductor 
as shown in Fig. [l|. The distance between dots is Tq. 




To 



FIG. 1. The triangular lattice with lattice constant To. 
Ti and T2 are the primitive lattice vectors. Ti = Tqx, 
T2 = \Tqx + ^Toy. The asterisks indicate the centers of 
the triangular plaquettes defined by the dots. 



Since the London theory cannot be solved in closed 
form for our system, we discretize the space of the unit 
cell including the region of the dot by dividing each edge 
of the unit cell into N intervals. The locations of the 
N x N grid points are given by 

X nu n, = ^fi + ^f 2l 0<ni J n 2 <JV. (1) 

The corresponding reciprocal lattice vectors are 



G, 



mxWi + m 2 W 2 , < mi,m 2 < N, (2) 



where W\ and W2 are the primitive vectors of the re- 
ciprocal lattice. 

Having discretized the space and defined the grid 
points, we now construct the current network of the unit 
cell. First, we define 

• Nodes: located at the center of each triangular 
plaquette formed by the three closest grid points. 
These nodes form a honeycomb lattice. 

• Links: are the edges of the honeycomb lattice. The 
link midpoints are located at the intersection of the 
links and the plaquette edges. 



We model our system as a wire network. Hence, the 
currents J S (X) only flow on the links. [See Fig. || and 
Fig. §.] 




FIG. 2. The links form a honeycomb lattice, and the cur- 
rent only flows on links. 



Fig. U illustrates a detailed section of our wire network. 
Three current links are associated with each grid point 
X, which may be written as 



M x ) 



Jsl(X) 

J s2 (X) 

J s3 (x) 



(3) 



where J s i(X) is centered at X+^ai, J S 2(X) at X+^a 2 
and J S 3(X) at X + i(<?i + a 2 ), with ai = j^Ti and a 2 = 



-k T 2 
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FIG. 3. Currents are depicted for the point X. Current 
conservation requires that J s i = J S 2 + J S 3- 
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We now enforce current conservation at the nodes. 
This gives two equations per grid cell (i.e., the cell formed 
by a! and 02). 



and 



J sl (X) = J S 2{X) + J s3 (X) 



J s3 (X) + J s2 (X + ai) = J sl (X + a 2 ). 



(4) 



(5) 



After Fourier transforming Eq. (0) and (^|) , we find that 
the current can be written in Fourier space as 



Js(G) 



I _ e -iG-S 2 
g—iG-a? g— iG-ai 



<f(G), G^O. (6) 



Since there is no net current flowing in the supercon- 
ductor, we have to include the G = component. There- 
fore, the current at the grid point X can be written as 
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We can think of 



Joi 
J02 
$(Gi) 
*(G 2 ) 

^(Gtv-i) 

as a vector, whose entries are determined by the cur- 
rent distribution in the unit cell. Since we are considering 
systems with antidots, one must have no current flowing 
in the defects. As will be discussed in greater detail in the 
following section, we impose this no-current constraint 
by expressing the current as a linear combination of ele- 
ments of a specially constructed basis. The elements of 
this basis are themselves current distributions satisfying 
the no-current constraint and are conveniently expressed 
via the vector formalism mentioned above. The basis ele- 
ments span the set of all current configurations for which 
current is excluded from the defects. 



III. THE METHOD 

Generally speaking, one calculates current distribu- 
tions by minimizing the London free energy subject to 



certain constraints, including the absence of current in 
the defects and a fixed number of vortices per unit cell. 
This shall be our program as well. Therefore, we first 
write the free energy per unit cell of our superconduct- 
ing system in the London approximation, 
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Y,\W)\ 2 + ^- c Y,MX)-A tot (X). (8) 
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In a thin-film superconductor, if the thickness is much 
smaller than the London penetration depth (i.e., d << 
A), the volume current density J is nearly constant in the 
thickness, so it is appropriate to write J s = Jd and de- 
fine X e ff — X 2 /d. The first term in Eq.(||) is the kinetic 
energy of the supercurrent, and the second term is the 
magnetic field energy, which is the sum of the contribu- 
tions of the supercurrent itself (J S (X) ■ A S (X)) and the 
interaction between the supercurrent and the external 
field (J S (X) ■ A ext {X)). However, symmetry arguments 
indicate that the latter term vanishes in the case of thin- 
film superconductors. Since the currents are localized 
at the link midpoints, we only need to evaluate the dot 
product in Eq.(||) at the three midpoints in each grid cell. 
Hence, the free energy can be written as 
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(9) 



In the rest of this section, we will briefly outline our 
method for solving the London theory for our system. 



A. Constructing the current basis 

To generate a basis satisfying the no-current condition 
in the defect region, we consider a functional which pe- 
nalizes currents in this region. A simple choice is 



W[J S ] 



v y,\j s (x )\\ 



1 / 

x„ 



(10) 



where the X^s are the points inside the dot, and Vo 
is a constant. We want to minimize W[J S ] subject to 
the constraint that the current in the unit cell is nor- 
malized. We can implement this constraint by adding 
a Lagrangian multiplier to H^Js]. The result is a new 
functional C, 
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£ = V ^|J s (X )| 2 -A{^|J s (X)| i 

X X 



1}. (11) 



If we write the supercurrent as a Fourier sum, plug 
Eq.(^J) into Eq.jll]), and carry out the sum in real space, 
we obtain £ is a function of joi,jo2 and ^(G). Reality 
of J S (X) implies that $(— G) = <f>*(G). This symmetry 
reduces sums over G to sums over G > 0. Next, we write 
3>(G) = $ re (G) + i<5j m (G) and minimize C with respect 
to Joij 302, $re(G) and $i m (G). This gives us four linear 
equations, which can be written as a matrix equation as 
Eq.@ to solve for j 01 , j 02 , $ re (G) and $ lm (G). 



Joi 

302 
*re(#l) 



M[G, G'] 



(12) 



( J > re(G'(Ar2_ 1) / 2 ) 
*im(G'l) 

< fim(G'( A r2_ 1 )/ 2 ) 

ioi 

$re(Gi) 

< J 5 re(G( A r2_ 1 )/ 2 ) 
$im(G*i) 



where M[G, G'\ is an (N 2 + 1) x (iV 2 + 1) matrix, and 
e x A/Vo. 

After solving the eigenequation, we obtain (A" 2 + 1) 
eigenvalues. The eigenvalues exhibit a clear hierarchy. 
Fig. |] is an example for AT = 11 with R = hhnm. q 
values are very high relative to the other A^ 2 + 1 — q 
eigenvalues. The large eigenvalues correspond to eigen- 
states for which currents within the defect are nonzero. 
Excluding these eigenstates allows us to construct a ba- 
sis satisfying the no-current condition. This can be done 
by inserting the entries of the eigenvector into Eq.(^)., 
yielding 3 x (A^ 2 + 1 — q) real space versions of the ba- 
sis elements J n (X) for n = 1,2, (N 2 + 1 — q), and 
j = 1,2,3. Any component of the current distribution 
can be expressed as a linear combination of the basis el- 
ements, i.e., J 3 (X) = J2n v nJi( x )- 



FIG. 4. The 122 eigenvalues for N = 11. There are six 
high values, corresponding to eigenvectors for which current 
is nonzero in the defect region. To = 600nm, R = 55nm. 



B. Calculating the vector potential 

Having constructed the current basis, we can calcu- 
late the vector potential associated with each element of 
the current basis. In the discretized scheme, the vector 
potential is 



allspace 
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(13) 



In evaluating Eq.([13|), one must include the contribu- 
tions from all three types of current. Using the fact that 
J(X + di) — J(X), where 3? is a lattice translation vector, 
we can write 
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ujj is the displacement from the grid point X to link 
j = 1,2,3 inside the grid cell. To evaluate the slowly 
converging sum in Eq.(nJ), we use the Ewald trickta. 



Defining £ = | X - X' + ujj - w- 

1 2 f°° , f 2,2 

£ Jo 



die 



- SR| , one uses the identity 



-4= dte^ 2 



-er/(£>7) + jerfc(£_r]) 



(15) 



to break up the sum into a long-ranged piece (the first 
integral) and a short-ranged part (the second integral). 
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The short-ranged term falls off rapidly for £ >> 1/rj, and 
is readily computed. On the other hand, the long-ranged 
term converges slowly in real space. However, with the 
aid of two-dimensional f pr . m .of Ewald's generalized theta 
function transformationML-J 



d2 g e - tl x+a 



-c^.,1 2 i(G+K)-X 



t 



-HG+iq-icSj-cs,,) 



(16) 



the long-ranged term can be recast in a form that 
falls off rapidly with increasing K in reciprocal space, 
making the sum tractable. (In this context, K = 
N(u\Wi + M2W2), where u\ and u 2 are arbitrary inte- 
gers.) The expression for A°(X), formed by combining 
the short and long ranged pieces of the sum, converges 
reasonably quickly for an optimal choice of 77. The to- 
tal vector potential in the system due to supercurrents 
at link j of X is the linear combination of A J n (X), i.e., 
A j (X) = Y,n v nAi{X), where A{(X) is the vector po- 
tential associated with the nth element of the current 
basis j£(X) via Eq.(pl). 



C. Minimizing the free energy in the presence of 
constraints 

The Ginzburg-Landau theory yields a set of conditions 
that allow us to compute the currents and fields. We 
will see that these quantities minimize the London free 
energy while satisfying system constraints. In the limit 
that London theory is quantitatively valid, the current 
has the form as 



— * p n — * p —* 

2m c 



(17) 



where e* is the superconducting electron charge, m* 
the mass, and n s the number of superconducting elec- 
trons per cm 3 . <p(r) is the phase of the order parameter, 
and is therefore a single- valued function. Singularities in 
f(r) at fi correspond to vortices: 



(18) 



where pi is the winding number for a vortex at fi 
Using the expression for J if) given in Eq.(|l7j), we find 



Ar) ■ dl= E -^{h J> V<p(f) • dl- -M 
2m* J c 
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where 4>o(= ch/e*) is superconducting magnetic flux 
quantum, and $b the magnetic flux contained in path. P 
is an integer which may arise in <f V(p(r)-dl for geometries 



in which the superconductor is not simply connected. For 
example, integrating around a dot can give § V(p(f) - dl ^ 
without any singularities. Minimizing the London free 
energy yields 



h(f) + A 2 V x V x h(r) 

h(r) + -^V x J(f) = 0. 
c 



(20) 



Hence, for any path that does not enclose a net vortex 
charge, we have 



J(r) • dl 



4ttA 2 



A(r) ■ dl, 



(21) 



precisely as given by Eq.(|l9|). Thus, it is sufficient 
to seek currents satisfying Eq.(|l9|). By symmetry, the 
integral § cell J(f) ■ dl vanishes if the integration contour is 
taken to be the unit cell boundary. From this condition, 
we see that 



(22) 



where <& uc is the flux through the unit cell. Eq.(|22 
tells us that the magnetic field can only have the same pe- 
riodicity as the dots if the number of flux quanta through 
a unit cell is an integer. For a thin-film superconductor, 
it is useful to write Eq.(19) in terms of the surface current 
density J s and ^ e ff- 



J S (X) ■ dl= h f - + P}. (23) 



Eq.(E3|) gives us two constraints to enforce: 



and 



UX) ■ dl= Y^r—iP ~ (24) 
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Within our network model, § J S (X) ■ dl takes the form 
of a discrete sum which, when evaluated, depends only 
on the basis expansion coefficients v n . The path of inte- 
gration is the set of links making up the boundary of the 
defect, and § Dot J s (X) ■ dl is the sum of the inner product 

of the currents J| on the links and the links ^o- I n the 
case of Eq. (|25| ) , the integral about a singularity, the path 
consists of the six edges of the hexagon centered at X, 
and § Sing J S {X) ■ dl is the sum of the six inner products 

J^lo- Fig. ^depicts an example of the path of integration 
in the discrete scheme. In the figure, the path of integra- 
tion about the defect contains 18 links, while paths about 
singularities are hexagonal, containing 6 links. 
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FIG. 5. The line integral paths around the dot located at 
C and a singularity located at Xi. Paths live on links. The 
border of the path around the dot consists of the 18 links 
making up the border of the defect. The path around the 
singularity is simply the border of the hexagonal plaquette 
centered at the singularity. 



the magnetic flux enclosed by the path of integra- 
tion, is 



$ B = / d 2 X ■ (h ext (X) + h s (X)) 

d 2 X-h ext (X) + <fdl-A s (X). 



(26) 



The vector potential A S (X) is calculated via Ewald 
method and is written in terms of J S (X). In a manner 
similar to that described above for the current line in- 
tegral, one can compute the line integral § dl ■ A S (X). 
Again, the result depends only on the coefficients v n . It 
is easy to calculate the flux due to the external magnetic 
field. One readily solves Eq.(22) for the applied field h ext , 
finding 



— (2^ 



p) 



(27) 



where a uc is the area of the unit cell. Hence, we have 
for the external magnetic flux, 



d 2 X-h ext {X) 



Qpath 



> £> + P), 



(28) 



where a pa th is the area enclosed by the line of integra- 
tion. Inserting the right side of Eq.(28) into Eq.(|24|) and 
performing the discrete version of the line integral yields 

N 2 +l-q 3 

kc 2 



47te*A 



eff 



-ir-N EiP ^ P ). (29) 



C locates the center of the defect. 7 J (G) is the sum of 
Jl ■ Iq, and varies with the paths. Nq is the number of 
plaquettes enclosed by the path around the dot. 

Similarly, the discrete version of Eq. (p5| ) is 
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where 77 1 = 1 



if 
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and rf 
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Note that Eq.(30) holds for each of the plaquettes ex- 
ternal to the dot, yielding a total of N 2 — (q+1) equations. 
Equations(^) and ( pO] ) allow us to solve uniquely for the 
coefficients v ni permitting one to calculate the currents 
corresponding to a particular vortex configuration. 



D. Solving for v n 's 

Since the current in the system is given by J° S {X) = 
J2n v nJi{X), solving Equations (p9|) and ( |30| ) for the co- 
efficients v n is tantamount to finding the current distri- 
bution. As previously mentioned, we seek N 2 + 1 — q of 
the coefficients v n , where q+l is the number of plaquettes 
comprising the defect. Eq.(p9|), and Eq.(30) represent a 
total of N 2 — q equations. However, it may be shown 
that only TV 2 — (q + 1) of these are independent. 

Requiring that the net supercurrent vanish yields 



and 



J2 v "3o2 = 0. 



(31) 



(32) 



One now has the desired TV 2 + 1 — q independent linear 
equations for the coefficients v n . In matrix form, we have 
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The coefficients a 
ing Eq.Q yields the 
Eq.(0), the free energy. 



are found from Eqs.([29| 
current distribution 



2|). Solv- 
and, via 
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IV. RESULTS 

A. Pinning Potentials 

The freedom to control the locations of vortices in the 
unit cell permits the calculation of pinning potentials. 
The simplest case, in which a vortex interacts with a 
pinning center containing zero flux quanta, is displayed 
in Fig. ^. One can see that the vortex is strongly pinned 
by the dot because the potential drops discontinuously 
when the vortex enters the defect zone. The potential 
surface in the vicinity of the dot forms a potential well 
with a length scale set by the defect radius. 




is qualitatively very similar to that discussed in Ref. |T^. 
As the darkened regions in Fig. || indicate, the intersti- 
tial pinning potential is quite shallow. Therefore, one ex- 
pects a higher mobility for intcrstitially pinned vortices 
than for vortices captured by the dots. 
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FIG. 7. Pinning Potentials for a single vortex in the unit 
cell of antidots with three different radii. To = 600nm, 
X e f f — 90nm. 



FIG. 6. Energy surface for a single vortex in the unit cell 
of antidots. T = 600nm, R = 0.06T , A = 1.5T . 



Examining a cross section of the pinning potential al- 
lows one to infer the pinning energy, the energy required 
to dislodge a single vortex from the interior of the defect. 
Fig. fj] displays potential cross sections for defect radii 
spanning the range between the coherence length £(T) 
and X e f f . A monotonic lowering of the pinning potential 
with increasing dot radius is evident from the plot. The 
graph in Fig. ^ supports the conclusion that the pinning 
energy initially increases rather swiftly as the radius of 
the defect increases beyond the coherence length. 

Interesting physics arises when one computes the pin- 
ning potential for a vortex interacting with a defect con- 
taining a single flux quantum. The pinning potential for 
a typical case is shown in Fig. |8|. The energy barrier sur- 
rounding the dot is perhaps the most salient qualitatively 
new feature compared to Fig. [?]. The barrier tends to re- 
pel the external vortex from the defect zone. However 
if the vortex manages to cross the barrier, the pinning 
potential abruptly decreases and the vortex finds itself 
in a local minimum. 

As the image indicates, another local minimum occurs 
at sites located at the centers of the equilateral trian- 
gles formed by the antidots. This interstitial minimum 
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FIG. 8. Energy surface in the unit cell of saturated antidots 
(n s = 1). To = 600nm, R = 0.067b, A = 1.5T . 



Fig. ^ displays a cross section of the single quantum 
pinning potential. The cross section is calculated along 
line A as illustrated in Fig. [I]. The domain in the graph 
spans the distance between adjacent defects. Evidently, 
the presence of an energy barrier surrounding the dot ne- 
cessitates the existence of an interstitial minimum. The 
inset offers a comparison with the case in which the pin- 
ning centers are not occupied by vortices. The pinning 
potential cross sections are qualitatively very similar to 
potentials calculated in Ref. and Ref. [n]. 
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FIG. 9. Energy of vortex vs. the distance from the center 
of the antidot along line A in Fig. |l|when each dot is saturated 
with one vortex. The inset shows both curves for empty and 
saturated dots. 



1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

R/R 1 



FIG. 10. Pinning phase diagram with weak dependence 
on A e //. The close-face squares represent the case of 
one captured vortex and one interstitial vortex, while the 
open-face triangles indicate that two vortices are captured. 
T = 600ran, 7?i = 0.05T . 



Fig. |l0| displays a pinning phase diagram. In the di- 
agram the saturation number n s , the number of flux 
quanta captured by the dot, is shown as a function of 
defect radius R and A e / / . The phase diagram is plotted 
for fixed £(T). Since the graininess of the discretization 
makes the notion of "radius" ambiguous, we take the de- 
fect radius to be the radius of the circle whose area is 
the same as the area enclosed by the defect boundary. 
A salient feature of the phase diagram is the absence of 
a strong dependence on the effective penetration depth, 
X e f j£9. In fact, the chief factor in the pinning phenomena 



for fixed £(T) seems to be the defect radius. The depen- 
dence of n s on R fits one's intuition. For small values 
of R, it is energetically favorable for the dots to capture 
one vortex each, while the remaining vortex resides at 
the interstitial site. For larger defects, the two vortex 
capture scenario becomes increasingly favorable. When 
the dot is made larger than about 0.15 lattice constant, 
the trapping of two quanta prevails and n s = 2. 



B. Pinning potentials for noncircular defects 

It is interesting to examine pinning phenomena for 
cases in which the pinning center has an elongated struc- 
ture. With the geometry of the triangular lattice, it is 
convenient to study defects in the shape of a rhombus, as 
illustrated in Fig. [El]. It is informative to calculate pin- 
ning potentials in which a vortex interacts with defects 
containing a single flux quantum. The potential surface 
is depicted in Fig. [l^. As for the circular case, energy 
minima appear near the centers of equilateral triangles 
formed by the defects. The structure of the interstitial 
sites, including the shallowness of the associated minima, 
is very similar to that of the circular case. As in the case 
of circular defects, the occupation of the rhombic defect 
by a vortex is associated with an energy barrier. How- 
ever, the barrier in the rhombic case is distinguished by 
a marked anisotropy. The barrier near the acute-angled 
corners of the rhombus is lower than the barrier height at 
any other part of the defect boundary. One readily con- 
cludes that vortices tend to enter or leave defects through 
the sharper points. This phenomenon is reminiscent of 
the manner in which electrons on a conductor preferen- 
tially discharge from the region of greatest curvature. 
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FIG. 11. The triangular lattice of antidots in the shape of 
rhombus. Each side of the dot has length S. 
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It is useful to compare the second capture energies of 
circular and rhombic defects of comparable area. (In this 
context, the TVth capture energy is defined to be the shift 
in energy which occurs when a vortex is moved from the 
interstitial region to the interior of the defect to join N— 1 
previously captured quanta.) To aid in the comparison, 
we display on the same plot pinning potential cross sec- 
tions for defects with both circular and rhombic shapes 
in Fig. [[3] and Fig. [l4| The lone circular defect pinning 
potential is plotted along a radial line for the circular 
case. The 3 rhombic potentials shown on the same graph 
correspond to paths of approach which the reader can 
identify from Fig. [fl| "Line 1" represents the easiest 
line of approach for the rhombus, in which one moves 
along the symmetry axis of the rhombus directly toward 
the sharpest corner of the defect. Two features of the 
potential are favorable for pinning. The "energy of ap- 
proach" , defined as the energy needed to move the vortex 
to the crest of the energy barrier, is lower for the rhom- 
bic defect than for the circular dot. The potential also 
indicates a greater pinning energy for rhombic defects. 
In fact, Fig. [l4| displays a case in which it is favorable 
for a rhombic defect to capture two flux quanta while 
its circular counterpart is only able to capture a single 
vortex. The lower energy of approach ensures faster ki- 
netics, hence swifter vortex capture for elongated pinning 
centers. The greater pinning energy, on the other hand, 
ensures that the vortices remain trapped once captured. 
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FIG. 12. Energy surface for a single vortex in the unit 
cell of saturated rhombic defects containing one vortex each 
(n s = 1). To = 600nm, S ~ 0.16T , A = 1.5T . 
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FIG. 13. Energy of vortex vs. the distance from the cen- 
ter of the antidot along line A in Fig. [l], line 1, line 2 and 
line 3 in Fig. |ll] when each dot is saturated with one vortex. 
To = 600nm, X eff = 1.5T , S = 2R ~ 0.16T . 
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FIG. 14. Energy of vortex vs. the distance from the cen- 
ter of the antidot along line A in Fig. [I], line 1, line 2 and 
line 3 in Fig. |ll] when each dot is saturated with one vortex. 
T = 600nm, A e// = 0.1T , S = 2R ~ 0.16T . 



V. SUMMARY 



In the work described in this report, we have used a 
wire network model to provide a detailed treatment of 
pinning phenomena in nanoscalc defect arrays. Our re- 
sults are consistent with experiments involving thin-film 
superconductors and are qualitatively very similar to the- 
oretical work on three dimensional superconductors con- 
taining columnar defects. 
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A comparison of pinning potentials suggests that elon- 
gated pinning centers have superior pinning characteris- 
tics to similar sized circular defects. 

Future work will extend this study to more general 
defect types and lattice geometries. We have fashioned 
techniques that permit the study of films with nontrivial 
three dimensional structures. Many of the experiments 
performed on nanoscale defect arrays are neither deep 
in the thin-film limit nor in the bulk limit. To model 
these systems properly, it is necessary to go beyond the 
thin-film approximation and provide a detailed three di- 
mensional treatment of the superconducting films. Such 
studies will be presented in future work. 
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